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Many human social phenomena, suchh as cooperation , the growth of settlements , 
traffic dynamics 0-0] and pedestrian movement |[7|-p!0| , appear to be accessible to mathemat- 



ical descriptions that invoke self-organization 1T1||T3 • Here we develop a model of pedestrian 
motion to explore the evolution of trails in urban green spaces such as parks. Our aim is 



to address such questions as what the topological structures of these trail systems are [|13 
and whether optimal path systems can be predicted for urban planning. We use an 'ac- 
tive walker' model |ll^-|T9[] that takes into account pedestrian motion and orientation and 
the concomitant feedbacks with the surrounding environment. Such models have previously 
been applied to the study of complex structure formation in physical |[T^-[16|, chemical 



and biological |1T8|,|19| systems. We find that our model is able to reporduce many of the 



observed large-scale spatial features of trail systems. 
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Previous studies have shown that various observed self-organization phenomena in pedes- 
trian crowds can be simulated very realistically. This includes the emergence of lanes of 
uniform walking direction and oscillatory changes of the passing direction at bottlenecks 
||7|JT0[] . Another interesting collective effect of pedestrian motion, which we have investigated 
very recently, is the formation of trail systems in green areas. In many cases, the pedestri- 
ans' desire to take the shortest way and the specific properties of the terrain are insufficient 
for an explanation of the trail characteristics. It is essential to include the effect of human 
orientation. To simulate the typical features of trail systems, we have extended the afore 
mentioned model of pedestrian motion to an active walker model by introducing equations 
for environmental changes and their impact on the chosen walking direction. 

First, we represent the ground structure at place r and time t by a function G{f,t) 
which reflects the comfort of walking. Trails are characterized by particularly large values 
of G. On the one hand, at their positions r = ra(t), all pedestrians a leave footprints 
on the ground (e.g. by trampling down some vegetation). Their intensity is assumed to 
be /(r)[l — G{f,t)/Gj^s,^{r)], since the clarity of a trail is limited to a maximum value 
Graax{r)- This causes a saturation effect [1 — G{r,t)/Gmax{r)] of the ground's alteration by 
new footprints. On the other hand, the ground structure changes due to the vegetation's 
regeneration. This will lead to a restoration of the natural ground conditions G'o(^) with a 
certain weathering rate 1/T(r) which is related to the durability T(r) of trails. Thus, the 
equation of environmental changes reads 

j:6{f-ut)), (1) 

a 

where S{r — fa) denotes Dirac's delta function (which yields only a contribution for f = fa)- 
The attractiveness of a trail segment at place r from the perspective of place fa decreases 
with its distance ||^— "^0(^)11 and depends on the visibility cr{fa)- Considering this by a factor 
exp(— ||r — fa\\/cf{fa)) and taking the spatial average by integration of the weighted ground 
structure over the green area, we obtain 

Vtr{fa,t) = J c/Ve-ll"^-^^"ll/"("^")G(f,t). (2) 
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The trail potential Vtr{ra,t) reflects the attractiveness of walking at place fa- It describes 
indirect long-range interactions via environmental changes, which are essential for the char- 
acteristics of the evolving patterns . 

On a plain, homogeneous ground, the walking direction of pedestrian a is determined 
by the direction of the next destination d^, i.e. ea{ra) = {da — ra)/\\da — ra\\- Without a 
destination, a pedestrian is expected to move into the direction of the largest increase of 
ground attraction, which is given by the (normalized) gradient Vf^Vtr{ra,t) of the trail po- 
tential. However, since the choice of the walking direction is influenced by the destination 
and existing trails at the same time, the orientation relation 

^ ,x da — + \^f^VtT:{ra,t) . . 

e-aVTait) = ^ (3) 

Wda-ra + V^M^aM 

was taken as the arithmetic average of both effects. Considering cases of rare interactions, 
the approximate equation of motion of a pedestrian a with desired velocity w° is 

= Vleai^a^t) . (4) 

A comparison of simulation results with photographs shows that the above described 
model is in good agreement with empirical obvervations. In particular, the evolution of the 
unexpected 'island' in the middle of the trail system in Figure ^ can be correctly described 
(Figure The goodness of fit of the model is quite surprising, since it contains only two 
independent parameters k = JT/o"^ and A = V^T/a, where denotes the average of the 
desired velocities This can be shown by scaling the model to dimensionless equations. 
The parameter A was kept constant. 

Our simulations base on a discretization of the considered area in small quadratic el- 
ements of equal size, which converts the integral (0) into a sum. Temporal and spatial 
derivatives are approximated by difference quotients. The presented examples begin with 
plain, homogeneous ground. All pedestrians have their own destinations and entry points, 
from which they start at a randomly chosen point in time. In Figure (Figure 0) pedestrians 
move between all possible pairs of three (four) fixed places. While in Figure ^ the entry 
points and destinations are distributed over the small ends of the ground. 
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At the beginning, pedestrians take the direct ways to their respective destinations. How- 
ever, after some time they begin to use aheady existing trails, since this is more comfortable 
than to clear new ways. By this, a kind of selection process pO|j2l|JT6|| between trails sets 



in: Frequently used trails are more attractive than others. For this reason they are chosen 
very often, and the resulting reinforcement makes them even more attractive. However, 
the weathering effect destroys rarely used trails and limits the maximum length of the way 
system which can be supported by a certain rate of trail usage. As a consequence, the 
trails begin to bundle, especially where different trails meet or intersect. This explains, why 
pedestrians with different destinations use and produce common parts of the trail system 
(Figures || and |). 

A direct way system (which provides the shortest connections, but covers a lot of space) 
only develops if all ways are almost equally comfortable. If the advantage k of using existing 
trails is large, the final trail system is a minimal way system (which is the shortest way system 
that connects all entry points and destinations). For realistic values of k, the evolution of the 
trail system stops before this state is reached (Figure ^. Thus, k is related to the average 
relative detour of the walkers. We conjecture that the resulting way system is the shortest 
one which is compatible with a certain accepted relative detour. In this sense, it yields an 
optimal compromise between convencience and shortness. 

Therefore, we suggest to use the above model as a tool for urban planners and landscape 
gardeners, who have the dilemma to build most comfortable way systems at minimal con- 
struction costs. For planning purposes one needs to know the entry points and destinations 
within the considered area and the rates of usage of their connections. If necessary, these can 
be estimated by trip chaining models , which are also needed in cases of complex lines of 



access and sight. The effects of the physical terrain and already existing ways can be taken 
into account by the function Go{r). By varying the model parameter k, the overall length 
of the resulting trail system can be influenced (Figures ^ and In the same way, one can 
check its structural stability. Presently, we are evaluating typical parameter values of A and 
K, by comparison of simulation results with real pedestrian flows which are reconstructed 
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from video films by image processing. These values shall be used for designing convenient 
way systems in residential areas, parks, and recreation areas by means of computer simu- 
lations (Figure 0). We expect that such way systems will actually be accepted, since they 
take into account the route choice habits of pedestrians. 

In summary, the presented active walker model is able to describe the self-organization 
and the typical structural properties of human trail systems. It will be interesting to relate 
our model to the work on space syntax [^. Repulsive interactions between pedestrians can 



be taken into account by generalizing equation @ in accordance with the social force model 
of pedestrian motion [|7|,p!0|l. However, these are only relevant in cases of frequent pedestrian 
interactions, in which they lead to broader trails. 
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FIGURES 





FIG. 1. Between the straight, paved ways on the university campus in Stuttgart- Vaihingen a 
trail system has evolved (center of the picture). Two types of nodes are observed: Intersections 
of two trails running in a straight line and junctions of two trails which smoothly merge into one 
trail. 




FIG. 2. The structure of the emerging trail system (light grey) essentially depends on the 
attractiveness parameter k. If k is small, a direct way system develops (left), if k is large, a minimal 
way system is formed, otherwise a compromise between both extremes will result (middle) which 
looks similar to the trail system in the center of Figure ^. 
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FIG. 3. The graphics illustrate the places and walking directions of pedestrians by arrows. The 
trail potential Vtr {r,t) is represented by a color scale (red = small, blue = large values). Starting 
with a plain ground, the structure of the trail system changes considerably during the simulation. 
Initially, pedestrians use the direct ways (left). Since frequently used trails become more com- 
fortable, a bundling of trails sets in which reduces the overall length of the trail system (right). 
The resulting way system (whose asymmetry is caused by differences in the frequency of trail 
usage) could serve as a planning guideline. It provides a suitable compromise between minimal 
construction costs and maximal comfort. Moreover, it balances the relative detours of all walkers. 
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FIG. 4. 
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I used for walking (irHontrast to paved areas). The similarity between the simulation 



result (left) Hid the trail^ystem oHthe university campus of Brasilia (right, reproduction by kind 
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